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Abstract: The deformation modulus of rock mass can either 
be measured using in-situ field tests or estimated empirically. 
However, the empirical method is an attractive and open 
choice for researchers to estimate the deformation modulus of 
rock mass economically, efficiently and accurately and many 
researchers have proposed empirical equations for estimation 
the deformation modulus of rock mass based rock mass 
classification systems. In this paper attempt has been made to 
collect and review some these published empirical equations 
developed so far based on Rock Mass Rating, Geological 
Strength Index, and Q rock mass classification systems. These 
empirical equations based on each classification system are 
grouped on the based on similar parameters. All the equations 
within certain group for all groups were plotted to obtain a 
viable range for which the equation is assumed to be reliable. 
An artificial data is generated for each group separately to 
from valid data points of existing equations. Keeping in view 
the trend of data points in the scatter plot of individual group, 
Gaussian and Sigmoidal type mathematical functions are 
fitted to obtain generalized equation for each group that 
consisted two or more than two valid equations. The new 
equations fitted to the data are optimized by minimizing 
residuals using Microsoft Excel built in Add- in "Solver". It is 
observed that equations using sigmoidal function predict more 
precisely than equations using Gaussian function. 

Key Words: Deformation modulus; Rockmass; Estimation 
methods; Empirical; Gaussian; Sigmoidal; 
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Introduction 



The deformation modulus of rock mass is an important 
parameter and used as essential input used in many numerical 
tools codes for response of ground to stress after excavation of 
tunnels within rock mass. The parameter can be directly 
determined directly using different in situ test procedures 
however; it requires significant cost and involves tedious 
operation process. Furthermore the intrinsic errors and other 
associated problems provide a base to search for alternative 
methods to estimate the deformation modulus of rock mass 
indirectly and is an interesting issue for research in the field of 
rock engineering. Historically the deformation modulus of rock 
mass is estimated using either empirical method utilizing rock 
mass classification systems, rock quality and seismic velocity in 
rock mass or utilizing certain models as by Duncan and 
Goodman, 1968 and Kulhawy, 1978 (Benson, 1986). 
Deere, 1969 made the first attempt to estimate the deformation 
modulus of rock mass establishing correlation between modulus 
obtained from in situ tests and Rock Quality Designation 
(RQD). The correlation was poor as there was great scatter in 



the data. Continuing the line Bieniawski, 1978 felt that the 
intact rock properties influenced the rock mass properties and 
developed an equation based RQD and modulus reduction 
factor which is the ratio between in situ modulus and laboratory 
modulus of rocks. Recently some new correlations have been 
developed based on either modulus ratio and RQD (Zhang & 
Einstein , 2004) and (Jing, 2003) correlating intact rock 
modulus E( 9 intact rock strength a c , weathering degree WD of 
rock mass with RQD (Kayabasi, Gokceoglu, & Ercanoglu, 

2003) and (Gokceoglu, Sonmez, & Kayabasi , 2003). Some of 
equations based on such parameters are listed in Table 1. 
However, all these equations are not considered for analysis as 
these equations are based on the parameters that are integral 
component of the major classification systems. 

The empirical estimation of deformation of rock mass using 
rock mass classification systems is an interesting and open topic 
for many researchers and further development is in progress in 
this area. For the last three decades a great number of empirical 
equations were proposed for estimation of deformation modulus 
of rock mass using numerous rock mass classification systems. 
Bieniawski, 1978 proposed the first empirical equation to 
estimate deformation modulus of rock mass using in situ data 
and RMR classification system. The journey of development of 
empirical equations based on rock mass classification systems 
is still continued and large number of such equations based on 
RMR (Bieniawski, 1978), (Serafim & Pereira, 1983), 
(Nicholson & Bieniawski, 1990), (Mehrotra, 1992), (Kim, 
1993), (Mitri, Edrissi, & Henning, 1994) (Read , Richards, & 
Perrin, 1999), (Mohammad , 1998), (Ramamurthy, 2001), 
(Ramamurthy, 2004), (Sonmez, Gokceoglu, Nefeslioglu , & 
Kayabasi, 2006), (Chun, Lee, Seo, & Lim, 2006)(Mohammadi 
& Rahmannejad, 2010) and (Shen , Karakus, & Xu, 2012) , GSI 
(Hoek & Brown, 1997), (Hoek, Carranza-Torres, & Corkum, 
2002), (Gokceoglu, Sonmez, & Kayabasi , 2003), (Carvalho, 

2004) , (Sonmez, Gokceoglu, & Ulusay, 2004), (Galera, 
Alvarez, & Bieniawski, 2005), (Hoek & Diederich, 2006), 
(Beiki, Bashari, & Majdi, 2010) and (Ghamgosar, Fahimifar, & 
Rasouli, 2010) and Q system (Barton, 1983) (Grimstad & 
Barton, 1993), (Diederichs & Kaiser, 1999), (Ramamurthy, 
2001), (Barton , 2002) and (Ramamurthy, 2004) have been 
developed so far while attempt have been made to include RMI 
(Palmstr0m & Singh, 2001) in the list. Majority of the 
equations are best fit on data for which 

the equations were developed. The behavior of rock mass varies 
and it is almost impossible to obtained such an equation that 
can be used globally however, it is possible to obtain more 
reliable general equation by including more/ new data and 
modify the relationship based on apparent trend of the data. 
Table 1: Equation based on intact rock properties, RQD and 
Weathering Degree 
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Eq. 
No. 



Equation 



Reference 



0.13 



0.00 



= EiW 



Ei(l + RQDI 100) 
WD 



(Eil oc)(\ + RQDI 100) 
WD 



0.0186fi<2£> -1.91 



1.5528 



Kayabasi et al., 
2003 



Gokceoglu et al, 
2003 



Zhang and 
Einstein, 2004 



2 Evaluation of the existing Empirical Equations 

In this research equations Based RMR, GSI and Q classification 
systems are evaluated as given in Table 2. These existing 
equations are grouped according to the classification system on 
which they based and divided into groups on basis of 
independent variables incorporated in the relationship. In this 
study the empirical equations based on RMR, GSI and Q 
systems are grouped into ten different groups subject to the 
input parameters used and each group is analyzed separately. 
Each individual equation in the groups is reviewed in detail to 
check their reliability, application and limitations. 

2.1 Generation of data points and Optimization of 
parameters 

An artificial data is generated for each group separately from 
valid data points of existing equations instead of using the 
actual data of the individual equation. For selected groups 
where generalization is possible and required,20 data points for 
each equation are generated. The equations are trimmed for 
their valid range and a scatter data is generated for each group. 
Keeping in view the trend of data points in the scatter plot of 
individual group, Gaussian and Sigmoidal type mathematical 
functions are fitted to obtain generalized equation for each 
group that two or more than two valid equations. The new 
equations fitted to the data are optimized by minimizing 
residuals using Microsoft Excel built in Add-in "Solver". 

2.2 Prediction performance and generalization 
2.2.1 RMR group I: 

The first ever equation for estimation of deformation modulus 
of rock mass based on RMR classification system developed by 
Bieniawski, 1978. However, the equation is only applicable for 
good quality rock masses with RMR values greater than 50. 
Keeping in view the limitation of equation 6, Serafim and 
Pereira, 1983 developed a new correlation based on data from 
Bieniawski, 1978 and data collected by the originators of 
equation 7, having RMR in the range from 25 to85. The 
equation overestimates the deformation modulus rock mass at 
for lower range i.e. RMR< 10 (Mohammad , 1998) and upper 
range i.e. RMR > 90 (Hoek & Diederich, Empirical estimation 
of rock mass modulus, 2006). Mehrotra, 1992 presented a new 
best fit RMR based equation similar to that of Serafim and 
Pereira equation on his own data. Kim, 1993 proposed an 
exponential equation however the equation extremely 
overestimates deformation modulus of rock mass for RMR >85. 
An attempt has been made by Mohammad, 1998 to reduce the 
value of estimated deformation modulus equal to zero for RMR 
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=0 by subtracting a factor of 0.562 GPa from the original 
Serafim, et al., 1983 equation but equation overestimate the 
modulus for upper range of RMR. Read , et al., 1999 
developed a third power equation by keeping E rm /Ej =1 for 

RMR =100. Chun et al, 2006 obtained an exponential 
relationship between RMR and deformation modulus for weak 
rock masses. Based on in-situ data from different parts of Iran, 
a polynomial equation was developed by Mohammadi and 
Rahmannejad, 2010 but due to poor asymptote the equation 
over estimates the deformation modulus RMR < 25 and 
overestimate for RMR approaching to 100. Shen , et al., 2012 
fitted a Gaussian function on data obtained from trend of data 
Bieniawski, 1978, Serafim and Pereira, 1983 and (Stephens & 
Banks, 1989). All the equations of the group were plotted as 
shown in Figure 1 to highlight the performance of each 
equation The group comprises of 9 equations however equation 
7, 8, 9, 11 and 13 are trimmed upto invalid range of RMR for 
which the equations were seems to overestimate the 
deformation modulus of rock mass. Keeping in view the trend 
of the data two different type of mathematical functions i.e. 
Gaussian and Sigmoidal are fitted to the data points as shown in 
Figure 2. The parameters of the models using these types 
functions optimized by minimizing the sum of error square i.e. 
residual using Solver of Microsoft Excel. The optimized 
general equations are given by equation (38) and (39) with 
residual of 3638.57 and 3612.4 respectively. It is observed that 
fitting of sigmoidal function bitterly described the correlation 
between RMR and deformation modulus of rock mass. 
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Figure 1: Comparative Plot of empirical Equations belong to 
RMR Group - 1. 
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Figure 2: Fitting of equation (38) and (39) on data generated 
from equations of Group - 1 
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2.2.2 RMR Group - II: 

Similar to Group - 1, the evaluation and refinement of equations 
in Group -II is continued by working on either new data or 
trying different mathematical functions on old data. Nicholson 
and Bieniawski, 1990 for the first time adopted the concept of 
including intact rock young modulus in relating deformation 
modulus of rock mass with empirical rock mass classification 
system RMR. Soon the idea is duplicated by Mitri et al, 1994, 
Ramamurthy, 2001, Ramamurthy, 2004, Galera et al, 2005, 
Sonmez et al, 2006 and recently by Shen et al, 2012. The idea 
to extend the expression of deformation modulus in this group 
was initially taken through correlations with observed 
deformations from case history studies and predicted 
deformations from finite element analyses. Findings of Sonmez 
et al, 2006 that deformation modulus of highly quality rock 
mass with weak intact rock pieces is highly controlled by intact 
rock properties has reinforced the idea to include the intact rock 
properties especially in equations while correlating deformation 
modulus of rock mass with rock mass classification system. 
The prediction performance of the empirical equation are 
checked by plotting and analyzing the equations listed under 
Group - II, it is observed that equation proposed by 
Ramamurthy, 2004 overestimates the value of E rm /Ej for lower 

range of RMR as illustrated by Figure 3. Gaussian and 
Sigmoidal type mathematical functions are fitted to the data 
points and optimized using Solver as shown in Figure 4. The 
optimized general equations are given by equation (40) and (41) 
with minimum residual of 2.434 and 2.383 respectively. 
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Figure 3: Comparative Plot of empirical Equations belong to 
RMR Group - II 
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Figure 4: Fitting of equation (40) and (41) on data generated 
from equations of Group - II 
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2.2.3 GSI Group - III: 

This group consists of only two published equations in which 
variation of deformation modulus is described as function of 
only GSI. The first equation was proposed by Gokceoglu, 2003 
followed by Ghamgosar et al, 2010 using exponential functions 
however, the later extremely overestimate the deformation 
modulus of rock mass for GSI > 80. No further analysis is 
carried out due the reason that a) the Ghamgosar et al, 2010 
equation is an extension of the Gokceoglu, 2003 by adding new 
data b) GSI only describe the quality of rock mass in terms of 
block size and properties of discontinuities and missing some 
important intact rock properties that may also responsible for 
controlling the deformation modulus. 

2.2.4 GSI Group - IV: 

The need for including an intact rock property as other 
controlling parameter in correlating deformation modulus of 
rock mass with GSI soon felt by the developer of the 
classification system while replacing RMR with GSI. The first 
equation of this group was proposed Hoek & Brown 1997 for 
rock masses having uniaxial compressive strength < 100 MPa. 
Beiki et al, 2010 developed and model based on genetic 
programming approach for estimation of rock mass using 
trigonometric function but model results negative value at both 
the extremes of GSI. Equations belong to this group are not 
used for further analysis as a) although the former equation is 
limited to weak rocks, it extremely overestimates the 
deformation modulus of rock mass due its poor asymptote and 
modified versions of the equation have been developed b) the 
second equation is based limited data with GSI range from 26 
to 82 majority occurring in range of 45 to 65 and gives negative 
deformation modulus values when extrapolated to extremes. 

2.2.5 GSI Group - V: 

The basic controlling parameters of deformation modulus of 
rock mass in the equations of this group are assumed as a c j, 

Damage Factor D and GSI. Only a single published equation is 
available in this group is proposed by Hoek et al., 2002 keeping 
in view that in weak rock masses the besides the other two 
parameters the deformation modulus of rock mass is controlled 
by the intact rock strength. However, the proposed equation is 
limited for weak rock mass and overestimates the rock mass at 
upper range of GSI. 

Further analysis of this group is not carried out due to a) there 
only a single equation in this group and no other data is 
available for comparison and b) the equation is limited to rock 
mass having a c [ <100 MPa. 

2.2.6 GSI Group - VI: 

Keeping in view the fact the deformation modulus of rock is not 
controlled by intact rock properties and rock mass quality but 
also influenced by method of excavation. For rock masses Hoek 
et al, 2002 founded that the deformation modulus is greatly 
influenced by the damage occurred due blasting and stress 
redistribution for excavation and established an equation for 
estimation of deformation modulus of rock mass from GSI and 
incorporate factor D. The equation is however is limited to rock 
mass with intact rock uniaxial compressive a c [ > 100 MPa. 

Based on in-situ data from China and Taiwan, Hoek and 
Diederichs, 2006 extended the previous equation to all rock 
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mass assumed to be homogenous and isotropic and established 
new equation based on GSI, and D. Similar to many other 
equations this equation has good prediction performance. 

The prediction performance of the empirical equation are 
checked by plotting and analyzing the equations listed under 
Group - VI, it is observed that equation proposed by Hoek et al, 
2002 overestimates the value of E rm for upper range of GSI as 

illustrated by Figure 5 . Gaussian and Sigmoidal mathematical 
functions are fitted to the upper and lower bound of the 
generated data and optimized using Solver as shown in Figure 6 
and 7. Four boundary equations (equation 42 to 45), two for 
each function type are produced with minimum residual of 
204.53 and 28.53 for upper bound and 2.19 and 0.0 
respectively. Finally combining the upper and lower bound 
equations of each function type, two general equations 
proposed for all values of D i.e. from 0 to 1. The optimized 
general equations are given by equation (46) and (47) 



respectively. 
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Figure 5: Comparative Plot of empirical Equations belong to 
RMR Group - VI 
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Figure 6: Fitting of Gaussian and Sigmoidal functions to lower 
bound data generated from equations of Group - VI 
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Figure 7: Fitting of Gaussian and Sigmoidal functions to upper 
bound data generated from equations of Group - VI 
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2.2.7 GSI Group - VII: 

The equations in this group are based on GSI, Ej and D. 

Carvalho, 2004 for the first time proposed 
equation estimation of E rm incorporating 

controlling parameters. Based on the assumption that when GSI 
= 100, the ratio of the modulus ratios of rock mass and intact 
rock is equal to unity, Sonmez et al, 2004 proposed a new 
prediction model for estimation E rm . However in both the 

equations the damage factors is not incorporated properly and 
estimate the same value of deformation modulus for all value of 
D i.e. 0 to lwhen GSI = 100. 

Based on in-situ data from China and Taiwan, Hoek and 
Diederichs, 2006 extended his own equation and incorporated 
Ej in equation earlier presented by the authors in the same 

paper. This equation is assumed to be more authentic as a) 
based on larger in-situ tests database b) variety of rock masses 
are incorporated and c) The mathematical function selected 
truly represents the trend of the scatter data. However the use of 
equation can lead the user to wrong end due to a) the intact rock 
samples for Ef are not always taken from behind the in situ test 

site and are not always true representative of the rock mass and 
b) although guidelines for selection of D are available (Hoek et 
al, 2002) but it is difficult for inexperienced engineer to obtain 
values for damage between the given values in the guideline. 

The prediction performances of the empirical equations are 
checked by plotting and analyzing the equations listed under 
Group - VII are plotted as illustrated by Figure 8. Artificial data 
points are generated to determine best fit generalized equation. 
Gaussian and Sigmoidal mathematical functions are fitted to the 
upper and lower bound of the generated data and optimized 
using Solver as illustrated by Figure 9 and 10 respectively. Four 
boundary equations (equation 48 to 51), two for each function 
type are produced with minimum residual of 0.03 and 0.02 for 
upper bound and 0.044 and 0.013 for lower bound respectively 
are obtained. Finally combining the upper and lower bound 
equations of each function type, two general equations 
proposed for all values of D i.e. from 0 to 1. The optimized 
general equations are given by equation 52 and 53 respectively. 
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Figure 8: Comparative Plot of empirical Equations belong to 
RMR Group - VII 
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Figure 9: Fitting of Gaussian and Sigmoidal functions to lower 
bound data generated from equations of Group - VII 
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Figure 10: Fitting of Gaussian and Sigmoidal functions to upper 
bound data generated from equations of Group - VII 

2.2.8 Q Group - VIII: 

Equations based on Q as proposed Barton, 1983 and Grimstad 
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and Barton, 1993 are for typical rock mass and typical range of 
Q value and estimate negative value of E rm for Q < 1. The 
equation proposed by Diederichs and Kaiser, 1999 is not based 
on actual Q system. The value of Q' parameter is obtained from 
RMR. The equation estimates a very high value of deformation 
modulus for low range of RMR and very low value for upper 
range of RMR. Further analysis of the equations in this group is 
not preceded due to a) the first and last equations poorly define 
the deformation modulus b) only the second equation seems to 
be valid but for Q > 1 and c) the originator of the first equation 
have proposed a modified version but more general equation by 
adding the possible influence of intact rock properties. 

2.2.9 Q Group - IX: 

Two equation belong to this group are collected from literature 
both developed by the same researcher. Ramamurthy, 2001 
proposed an equation assuming that the deformation modulus is 
dependent not only on Q index of rock mass and included E( 
while establishing equation for E rm . From the analysis of the 
equation it is revealed that for Q = 1000, the reduction factor 
i.e. E rm /E[ reach just to value of 0.75. Later on Ramamurthy, 

2004 modified their own equation by adding new data using the 
similar mathematical function. Although for Q = 1000, this 
equation do not achieve the theoretical E rm /E[ = 1 , the equation 

is considered to be more realistic in this group and no further 
analysis for development of general equation is carried. 

2.2.10 Q Group - X: 

Based on practical experience and case histories, Barton, 2002 
proposed a modified version equation by including o c { while 

correlating E m with Q index. The equation is suitably fit to data 

on which it was developed. In this group only one equation is 
available and no further generalization is required. 

3. Conclusion 

Analysis of this research shows that most of the equations for 
empirical estimation of deformation modulus of rock mass from 
rock mass classification systems are based on RMR. However, 
due to wide acceptance of GSI, the attraction to develop 
new/modified equations based on the classification system has 
been increased. In- spite of the popularity of Q system the 
progress and development is negligible. The reason might be 
the fact that the former two classification systems have broader 
applicability in rock engineering while the latter is limited to 
tunnel engineering. Majority of the equations either 
overestimate or underestimate the modulus at certain end and 
the equations need to be truncated. 

Group III, IV, V, VIII, IX and X are either limited to single 
equation or the validity of equations in the groups is limited and 
not suitable for generalization. Subject to the constraints new 
general equations are developed for Group I, II, VI, and VII 
using Gaussian and Sigmoidal type functions. Eight Simplified 
equations using Gaussian and Sigmoidal functions are 
proposed. The proposed equations in this research are best fit 
equations but selection must be made subject to availability of 
controlling parameters authentic data. It is concluded from the 
residuals the equations using sigmoidal function give better 
prediction than equations using Gaussian function for all groups 
i.e. Group I, II, VI and VII. 
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Table2: Empirical estimation of rock mass deformation modulus using rock mass classification systems 



System 


Group 


input 


Eq. 


Mathematical 


Empirical equation 


Reference 








No. 


function used 







RMR 



RMR 



Linear 



F rfTt = 2RMR - 100 



Bieniawski, 1978 



10 



Inverse log 
Inverse log 
Exponential 
Inverse log 



Serafim and Pereira, 1983 
Mehrotra, 1992 



E rm = 300^™* 10" 



E m = 10 : ^-^d _ 0 562 



11 



Third power ^=0.1 (RMR / 1 0] 3 



Kim, 1993 
Mohammad, 1998 
Read et al, 1999 



12 



Exponential 



Chun et al, 2006 



13 



Polynomial 



E rm = 0.0003RMR 3 - 0.0193RMR 2 + 
0.315RMR + 3.4064 



Mohammadi and 
Rahmannejad, 2010 



II 



f.::f.-uu\ 



14 



Gaussian 



F = 



110 e - " 



Shen et al, 2012 



RMR 

&Ei 



15 



Complex 



16 



Trigonometric 



17 



18 



Exponential 
Exponential 



E rrK = 0. OIF, (0.0028 RMR- -0.9*===3 



Nicholson and Bieniawski, 
1990 



r = 

Y~¥Ti. 



/tsRMR\\ 

0 ' 5 -^', 1 - cos (-T5rJJ 



h = 



E { 3 [ ±-* " 



27* _ pr „ - D.DD35 [5 ( LCD - MjO" 



Mitriet al, 1994 



Ramamurthy, 2001 
Ramamurthy, 2004 



GSI 
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IV 



19 



20 



Exponential 
Power 



F = 



f r^£-- : : 



Galera et al, 2005 



r = 



( r.-^--.::; - 1 1 - ,=:.-: ■=: 



Sonmez et al, 2006 



21 



f R M R - i i 

Gaussian JT r m = 1 . 1 4 Ii> ~ ■ * * 



Shen et al, 2012 
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Exponential 



E m =0.145e D " D&4G5J 



Gokceoglu, 2003 



23 



Exponential 



E rm = 0.0912e DDlE66&S[ 



Ghamgosar et al, 2010 
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and 
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Power 
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Hoek and Brown, 1997 
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Beiki et al, 2010 
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Hoek et al, 2002 
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VI 


GSI 
andD 


27 


Power 


o, > 100 MPa 

For 


Hoek et al, 2002 
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Sigmoidal 




Hoek and Diederichs, 2006 
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29 
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Carvalho, 2004 
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Sonmez et al, 2004 
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Hoek and Diederichs, 2006 
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32 


Logarithmic 




Barton, 1983 








33 


Logarithmic 


Ern, = 25]0gQ 


Grimstad and Barton, 1993 








34 




where 

Q"=10(RMR-^| 

21 


Diederichs and Kaiser, 1999 
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c nr_ — c i c 


Ramamurthy, 2001 
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